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According to the CMB observations, Mielczarek f [4p|p evaluated the reheating temperature, which 
could help to determine the history of the Universe. In this paper, we recalculate the reheating 
temperature using the new data from WMAP 7 observations. Based on that, we list the approximate 
solutions of relic gravitational waves (RGWs) for various frequency bands. With the combination 
of the quantum normalization of RGWs when they are produced and the CMB observations, we 
obtain the relation between the tensor-to-scalar ratio r and the inflation index f3 for a given scalar 
spectral index n s . As a comparison, the diagram r — j3 in the slow-roll inflation model is also given. 
Thus, the observational limits of r from CMB lead to the constraints on the value of /3. Then, we 
illustrate the energy density spectrum of RGWs with the quantum normalization for different values 
of r and the corresponding /3. For comparison, the energy density spectra of RGWs with parameters 
based on slow-roll inflation are also discussed. We find that the values of n s affect the spectra of 
RGWs sensitively in the very high frequencies. Based on the current and planed gravitational wave 
detectors, we discuss the detectabilities of RGWs. 

PACS number: 04.30.-w, 98.80. Es, 98.80.Cq 



I. INTRODUCTION 



General relativity and quantum mechanics predict a stochastic background of relic gravitational 



waves (RGWs) [l|-|6j, generated during the early inflationary stage. The primordial amplitudes could 
be determined by the quantum normalization at the time of the wave modes crossing the horizon 
during the inflation. After that, the evolution of RGWs are mainly determined by a sequence of 
stages of cosmic expansion including the current acceleration 0, 0], since the interaction of RGWs 
with other cosmic components is typically very weak. Therefore, RGWs carry a unique information 
of the early Universe, and serve as a probe into the Universe much earlier than the recombination 
stage. As an interesting source for gravitational wave (GW) detectors, RGWs exist everywhere 
and anytime unlike GWs radiated by usual astrophysical process. Moreover, RGWs spread a very 
broad range of frequency, 10 -18 — 10 10 Hz, making themselves become one of the major scientific 
goals of various GW detectors with different response frequency bands. The current and planed 



GW detectors contain t he g round-based interferometers, such as LIGO [9|, Advanced LIGO 

day 



VIRGO 
10 2 - 10 3 



12, 13 , GEO 14 



15 , LCGT 16] and ET h 



AIGO 

Hz; the space interferometers, such as the future LISA 
frequency range 10~ 4 — 10 



Hz, BBO mj, 22 




aiming at the frequency range 
which is sensitive in the 



23 which both are sensitive in the 



, and DECIGO r 

frequency range 0.1 — 10 Hz; and the pulsar timing array, such as PPTA (24 , 25|| and the planned SKA 
~26| with the frequency window 10 -9 — 10~ 6 Hz. Besides, there some potential very-high-frequency 
GW detectors, such as the waveguide detector [27 



around GHz [28 
interferometers 



the proposed gaussian maser beam detector 
and the 100 MHz detector with a pair of 75-cm baseline synchronous recycling 
291 1 . Furthermore, the very low frequency portion of RGWs also contribute to the 



anisotropies and polarizations of cosmic microwave background (CMB) [3 Of, yielding a mag netic 
type polarization of CMB as a distinguished signal of RGWs. WMAP [31 
ground-based ACTPol j36j and the proposed CMBpol [37| are of this type. 



Planck 35], the 
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The reheating temperature, Trh, carries rich information of the early Universe, and relates to the 
decay rate of the inflation as Trh oc Recently, the reheating temperature was evaluated 



[40| according to the CMB observations by WMAP 7 in the frame of the slow-roll inflation [34 
Then, the expansion histories of different stages could be determined subsequently. In this paper, 
we reevaluate the reheating temperature using the latest observational data from CMB, and adopt 
the resulting expansion periods of different phases of the Universe as references. The referenced 
reheating temperature can help us to divide the phases of the Universe definitely. The evolutions 
of the RGWs at various phases can be determined subsequently, and the primordial amplitude 
was normalized due to the quantum condition during inflation [H, Hj|. For the present time, the 
solutions of RGWs can be obtained for different frequency bands corresponding to the modes re- 
entered the horizon at different phases. Note that, this model of RGWs is free from the slow- roll 
inflation. Therefore, the above reheating temperature based on the slow-roll inflation just serves as a 
reference. On the other hand, the anisotropics due to the tensor metric perturbations (gravitational 
waves) can be scaled to those due to the observations of the scalar perturbations by introducing a 
parameter r called tensor-to-scalar ratio. Combining the observations of the CMB and the quantum 
normalization of RGWs when they are generated, a constraint condition is arrived between the ratio 
r, the inflation index /3, and the index /3 S describing the expansion behavior of the Universe from 
the end of inflation to the reheating process. For the chaotic inflation with a quadratic potential 
V = ^m 2 4> 2 , which means j3 a = 1, the diagram r — j3 will be illustrated for given values of the scalar 
spectral index n s . The resulting spectra of RGWs for different values of r the corresponding f3 will 
be demonstrated. For comparison, we will also discuss the spectra of RGWs with the values of r 
and j3 predicted by the slow-roll inflation itself. To this end, the spectra of RGWs given by different 
models and different parameters will confront the various current and planed GW detectors. 

The outline of this paper is as follows. In Sec. II, we recalculate the reheating temperature using 
the latest data from CMB and plot it as a function of the scalar spectral index. Based on that, the 
scale factor a(r) is specified for consecutive stages of cosmic expansion. In section III, we present the 
resulting approximate solutions of the spectrum of RGWs for various frequency bands. In section IV, 
the spectra of RGWs for different values of parameters are shown and some comparisons between the 
spectra based on quantum normalization and those based on slow-roll inflation will be given. Some 
discussions are summarized in Sec. V. Throughout this paper, we use the units c = H — ks — 1- 
Indices A, /i, u,... run from to 3, and i, j, k,... run from 1 to 3. 



II. THE EXPANSION HISTORY OF THE UNIVERSE 

For a spatially flat (k = 0) universe the Robert son- Walker spacetime has a metric 

ds 2 = a 2 (r)[-rfr 2 + S^dxP], (1) 

where r is the conformal time, and the scale factor a(r) is described by the following successive 
stages 0,0: 
The inflationary stage: 

a(r) = 1 \t\ 1+i3 , -oo<T<7i, (2) 

where the inflation index (3 is an important model parameter. The special case of /3 = —2 corresponds 
the exact de Sitter expansion. However, both the model-predicted and the observed results indicate 
that the value of /3 can differ slightly from —2. 
The preheating stage : 

a(r) = a z \r - T p \ 1+f3 % n < r < t s , (3) 
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where the parameter f3 s is usually taken as a constant. After inflation, the inflation field undergoes 
coherent oscillations at the bottom of potential well. The reheating process takes place the Hubble 
parameter H falls to the value of the inflation decay rate T^, if we assume that the reheating is 
instantaneous. Then the preheating process lasts from the end of inflation to the happening of 
reheating. Below, we will discuss the value of (3 S in detail. 
The radiation-dominant stage : 



o(r) 



(t-t 6 ), t s <t<t 2 . 



The matter- dominant stage: 



o(r) = a m (r - r m ) 2 , r 2 < r < t e . 
The accelerating stage up to the present time To 0): 

o(r) = l H \r - rJ~ 7 , t e < r < r , 



(4) 
(5) 
(6) 



where 7 is a Q a -d ependent parameter, and Q\ is the energy density contrast. To be specific, we 
take 7 ~ 1.97 [4l| for = 0.73 34J in this paper. It is convenient to choose the normalization 
l r o — To I = I? i- e -5 the present scale factor a(ro) = Ih- From the definition of the Hubble constant, one 
has Ih = j/Hq, where Hq = 100 h km s _1 Mpc _1 is the present Hubble constant. We take h ~ 0.704 



34j throughout this paper. Imposing (3 and (3 S as model parameters, only 12 constants remain in all 



the expressions of the scale factor. The continuity of a(r) and a'(r) at the four given joining points 
T i, T s, T 2 an d r E provide 8 constraints. If we also know the expansion history of various stages, i.e., 
the definite values of Ci = a(j s ) / a{ji) , ( s = a(r 2 )/a(T s ), ( 2 = a ijE) / / o(r 2 ) , and ( E = a(T~o) / cl{t e ) , all 
the 12 constants can be fixed completely j42| . 

The late-time universe is well know. For the ACDM model, the density of dark energy, which 
drives the accelerating expansion of the universe, is constant. Based on that, one easily has ( E = 
1 + z E = (^A/^m.) 1 ^ 3 — 1-4, where z E is the redshift when the accelerating expansion begins. For 



the duration of the matter-dominated stage, we get straightforwardly (2 



with z, 



eq 



3240 34 



a(r ) a(r g ) 
a(r2) a(ro) 



'J + Zeq)Ci 



However, the histories of the radiation-dominated stage and the preheating 
stage are not known well. Recently, Mielczarek 40[ proposed a method to evaluate the reheating 
temperature, Trh, under the frame of the slow-roll inflation model combing the observations from 
WMAP. Using this method, one can easily obtain the following expression: 



7i 



15 mpi 



RH 



7/2 



T, 



CMB 



exp 



6 



1 



(7) 



where mpi = l/y/G = 1.22 • 10 GeV is the Plank mass, g+ s = 3.91 counts the effective number 
of photons plus three species of massless neutrinos contributing to the radiation entropy during the 
recombination [43[, n s is the scalar spectral index, A s = (2.43 ± 0.11) • 10~ 9 is the amplitude of 
the scalar power spectrum at the pivot physical wavenumber /cq = 0.002 Mpc -1 [44j, and Tcmb = 
2.725 K = 2.348 ■ 10~ 13 GeV is the present temperature of CMB. Note that, in deriving Eq. (ED, 
the approximation of (3 ~ — 2 was used, and a quadratic potential V((p) = \m 2 4> 2 of the scalar field, 
where the scalar mass m can be fixed to be ~ 10 13 GeV 



451 ] by CMB observations, was assumed, 
during the prehearting stage, when the oscillating frequency of is much 
greater than the expansion rate of the Universe, the coherent scalar field oscillations behave like a 
fluid with p = wp, where the equation of state w depends upon the form of the scalar potential V{4>). 
Say V((f)) = A0", one has w = ^| and p decreases as a ~ 6n ^ n+2 \ For V{<p) = ^m 2 (j) 2 , i.e., n = 2, one 



As analyzed by Turner [46 



has w = and p oc a . This result was also verified by Martin and Ringeval |47J using numerical 
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FIG. 1: The reheating temperature Trh as a function of the scalar spectral index n s , where p s = 1 and (?* s = 3.91 are set. 



method, and it was found that the average w never deviates from zero exceeding 8%. Hence, during 
the evolution of the prehearting stage, the energy density drops approximately in the same way as 

and the scale factor evolves as a(t) oc t 2 ^ 3 , which means (3 S = 1 
1 except that we write it explicitly. Moreover, the 



18 



the matter-dominated Universe 
from Eq. ([3]). So, in the following we set /3. 



following relations in the slow-roll inflation [49 



m 



pi 



16tt V V 



V 



m 



V 



pi 



">7T 



6e + 2rj, 



(8) 



were also used in deriving Eq. (J7J). According to the observations of WMAP 7, Mielczarek obtained 
T RH = 3.5 • 10 6 GeV with n s = 0.963 ± 0.012 and A s = 2.441+°;°| ■ 10~ 9 , assuming g M = 2 being the 
effective number of relativistic degree of freedom only contributed by pho tons during recombination. 
One can recalculate Trh using the updated data given by WMAP 7 [34J. For example, for WMAP 
Seven-year Mean, n s = 0.967 ± 0.014, one has Trh — 5.8 • 10 14 GeV with a relative uncertainty 
gt(Trh)/T rh « 77; and for WMAP Seven-year ML, n s = 0.966, one has T RH ~ 2.8 • 10 12 GeV. Thus, 
one can see that Trh depends on the value of n s very sensitively. This is because the expression of 
Trh contains a exponential factor exp (t^-), which is dependent on the value of n s very sensitively. 
For a more general potential with power-law form V = \<ft n , one can easily obtain the e-foldingor 
number A^bs — 2 (l-n ) us i n S Eq. (14) in Ref. jioj], and in turn the resultant Trh has an exponential 

factor exp ( 2{i-n\ ) • ^ or exam pl e ; n — 4, which implies w — 1/3 and (3 S = 0, Trh contains a factor 
exp (tj — ). Therefore, Trh also depends on the power index n, and in turn, on (3 S very sensitively. We 
plot the reheating temperature as a function of the scalar spectral index in Fig. [TJ From the bottom, 
the big bang nucleosynthesis (BBN) gives a constraint of the reheating temperature Trh ^ 10 MeV 
[Ho| . From the top, the constraint comes from the energy scale at the end of inflation Trh ^5 10 16 
GeV. Under these constraints, one can see that the scalar spectral index should locate at the range 
of 0.959 < n s < 0.967. 

After the instantaneous reheating, the universe is filled with the relativistic plasma. We assume 
the expansion of the relativistic gas to be adiabatic, which is valid until the entropy transfer between 
the radiation and other components can be neglected. Therefore, the conservation of the entropy 
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gives the increase of the scale factor from the reheating till the recombination 40 



rp / \ 1/3 

J- RH / 9*s 



a(r s ) T rec \g 



(9) 



where a rec and T rec stand for the scale factor and the temperature at the recombination, respectively. 
g* s counts the effective number of relativistic species contributing to the entropy during the reheating. 
In the standard model of elementary particles, one has g* = g^s = 106.75 at the energy scale above 
~ 1 TeV j43[ , where g* counts the effective number of relativistic species contributing to the energy 



density during the reheating. Moreover, as pointed in [40j, if the temperature of reheating is greater 
than the electroweak energy scale, Trh ~ 300 GeV, one may expect that g* > 106.75. Thus, in this 
paper we set g* = g* s = 106.75. Based on Eq. ([9]), one easily obtain 

^ = a(r 2 ) a rec = T RU ( g*s\ 1/3 ^ 



a-rec a[T s ) T CMB (1 + z eq ) \g. 

where we have used T rec = Tcmb{^- + z rec)- Finally, for the slow-roll inflation, the increase of the 
scale factor during the preheating stage with /3 S — 1 is given by (4o| 



15m 4 /L(l 



2\ V3 



64vr 2 ^T 4 H / 

Even though ( s and Ci are dependent on g*, T RH is independent on g m as shown in Eq. (JTj). For 
WMAP Seven-year Mean, n s = 0.967, one has ( s ~ 2.3 • 10 24 and Ci ^ 4.86; while for WMAP Seven- 
year ML, n s = 0.966, one has ( s ~ 1.11 • 10 22 and Ci — 6.07 ■ 10 3 . In the following, the resulting ( s 
and (i under the frame of the slow-roll inflation will serve as referenced tools to obtain the solutions 
of the RGWs, even though some calculations of RGWs based on quantum normalization are free 
from the slow-roll inflation. 

III. RGWS IN THE ACCELERATING UNIVERSE 

In the presence of the gravitational waves, the perturbed Robert son- Walker metric is given by 

ds 2 = a 2 (T)[-dr 2 + (5 i:i + hijjd.r'd.r 1 . (12) 

where the tensorial perturbation hij is traceless h % i = and transverse hijj = 0. It can be decomposed 
into the Fourier A;-modes and into the polarization states, denoted by a, as 

Mr,x)= £ fj^j-A^ir)^, (13) 



<T= + ,X 



where h^l* (r) = hfjf* (r) ensuring that be real, the comoving wave number k is related with the 
wave vector k by k = (Sijk 1 ^) 1 / 2 , and the polarization tensor ej^ satisfies 0]: 

^ = ^. 6^ = 0, #V = o, e l; ) (-k) = e ;; ) (k). (i4) 

In terms of the mode h!^\ the wave equation is 

ht )n {r) + 2?&hP'(T) + k 2 ht\r) = 0, (15) 



G 



where a prime means taking derivative with respect to r. The two polarizations of h k a \r) have 
the same statistical properties and give equal contributions to the unpolarized RGWs background, 
so the super index (a) can be dropped. Introducing a new notation [3, 0, 0]: Mas(t) = a(r)/i fc (r), 
Eq. (IT5"l) reduces to 



a" 



Vk+[k--)Vk = 0. (16) 

For the high-frequency limit, i.e., the term a" /a can be neglected, the solution to Eq. (I16|) has the 
usual oscillatory from: /ifc(r) = Akt~ %kr + Bke lkT , where the constants Ak and Bk are determined by 
the initial conditions, and hk decreases adiabatically with the expansion of the universe: 

±ikr 

h k « — -. (17) 
a(r) 

On the other hand, when the term a" /a is dominant, the dominant solution of Eq. (fl6|) is /ifc(T) = 
Cfca(r) for growing functions a(r) in expanding universes, and 

hk = const. (18) 

In other words, when the wavelength of the mode hk{r), A = 2ira/k, is much less than the horizon of 
the universe, 1/H = a 2 (t) / 'a' 'Or) , hk(r) decays with the expansion of the universe; while if A 3> 1/H, 
hk{r) will keep constant [l|, 0]. The wavelengths of various modes hk are stretched outside the 
horizon during the inflation stage, and they will be constant until re-enter the horizon again. When 
the mode hk re-enters the horizon, it will decay as hk oc l/a(r). 
The spectrum of RGWs h(k, r) is defined by 



00 2r ,dk 



(^(r,x)%(r,x)}= / h'(k,r)—, (19) 



o 



where the angle brackets mean ensemble average. The dimensionless spectrum h(k, r) relates to the 
mode hk(r) as 51 



h(k,r) = ^k 3 / 2 \h k (r)\. (20) 

7T 

The primordial spectrum of RGWs at the time Tj of the horizon-crossing during the inflation has a 
power-law form 0, 0, 0, : 



h(k,Tj=A^—) , (21) 

where the constant A representing the initial condition is to be determined both from theories 
and observations. An exact de Sitter expansion, i.e., (3 = —2, yields a scale-invariant spectrum. 
Consider Eqs. f lT7|) . fTT8|) and f l20|) . one knows that the spectrum h(k, r) will decay as oc l/a(r) when 
it reenter the cosmic horizon, 1/H. In the following, let us discuss the properties of the present 
RGWs spectrum h(k, tq) with different frequency bands. The characteristic comoving wave number 
at a certain conformal time r x is give by 



(22) 



FIG. 2: The evolution of a /a. 




(23) 



As shown in Figj2j for the comoving wave number k < ks, the modes of RGWs have been outside 
the horizon all the time. Thus, these modes never decay and keep their original amplitudes at 
present. For, fcg < k < kn, the modes entered the horizon before the beginning of the acceleration 
and went out the horizon again before the present time. We denote the time of the mode k entering 
and going out of the horizon as r* and r**, respectively. The corresponding scale factor are marked 
as a* and a**, respectively. Then, the present spectrum for fc^ < k < kn can be written as 



According to Eqs.([5]), ()6]) and (|22|) . one can easily obtain a^/aijE) = ikE/k) 2 and a{rE)/a** = 
(fc/Zc^) -7 , which lead to h(k,r ) = h(k,T i )(k/kE)~ ( ' 2+ ^ ■ For k > k H , after the modes entering the 
horizon, they have been inside the horizon up to now. Therefore, the spectrum for all the modes of 
k > ku has the following form: 



where we have still used the notes a* standing for the scale factor when the k— mode enters the 
horizon. Using the similar analysis given above, one can get the present spectrum h(k,To) for 
different bands of wave number which correspond to different stages of the universe. Generally, we 




(24) 




(25) 
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summarize the approximate solutions of h(k, Tq) uniformly as follows , 



h(k, t ) = A ( — ) , k<k E - (26) 



k 



(8-7 



2+7 



h(k,T ) = A[—j (l + z E )-—, k E <k<k H - (27) 
h(k 7 T )=A[—) (I + ze)-^, k H <k<k 2 ; (28) 



7 \ 1 + /3 /I \ 

= A(— J f-^J(l + ^)-^, k 2 <k<k s ; (29) 



h(k 1 T ) = A^—j ^j(i + ^)-— , k s <k<h. (30) 

2+7 

The factor (1 + z E ) 7 means a reduction due to the accelerating expansion of the universe. The 
above equations reduce to the corresponding results shown in where 7 = 1 was assumed. 

Another important quantity often used in constraining RGWs is its present energy density pa- 
rameter defined by flew =< Pg > / Pc where p g = 32^^.7,0^0 * s the energy density of RGWs, and 
p c = 3Hq/8ttG is the critical energy density. A direct calculation yields 0, [HJ 



fupper 
flou 



tt GW = I <W)y, ( 31 ) 



with 



^(/) = ^ c 2 (/)(^) (32) 

being the dimensionless energy density spectrum. We have used a new notation, h c (f) = h(f, r ) / y/2, 
called characteristic strain spectrum [5| or chirp amplitude jUJ . The lower and upper limit of integra- 
tion in Eq. (l3Tj) can be taken to be fi ow ~ f E and f upper — fi, respectively, since only the wavelength 
of the modes inside the horizon contribute to the total energy density. 
In the present universe, the physical frequency relates to a comoving wave number k as 

f = * = k (33) 

2na(r ) 2-kI h y ' 

Thus, from Eqs. (T231) and fl33|) one can easily get the each characteristic physical frequency: /# = 
Ho ~ 2.28 • 10~ 18 Hz, f E ~ 1.93 • 10" 18 Hz, f 2 ~ 9.3 ■ 10~ 17 Hz. Moreover, f s ~ 2.14 ■ 10 8 Hz and 
/1 ~ 4.71 ■ 10 8 Hz for n s = 0.967; f s ~ 1.03 • 10 6 Hz and h ~ 8.04 • 10 7 Hz for n s = 0.966. The 
values of f\ are below the constraint from the rate of the primordial nucleosynthesis, f\ < 3 x 10 10 
Hz When the acceleration epoch is considered, the constraint becomes /i~4x 10 10 Hz. Note 
that, Eq. ( l2"B"j) implies that fi depends on the value of fi s but not (5. Therefore, f\ is a fixed value 
once f3 s has been chosen. 



IV. OBSERVATIONAL CONSTRAINTS AND THE DETECTION 



As well known, the anisotropies and polarizations of CMB are contributed by two parts: tensor 
metric perturbations (gravitational waves) and curvature perturbations. Moreover, the B-mode 
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CMB polarization is only produced by tensor perturbations. In literatures [31n34(, the tensor-to- 
scalar ratio r is often introduced 

r = (34) 

where A^(/c ) and A^(fc ) are the power spectrum of the tensor perturbations and curvature per- 
turbations evaluated at the pivot wavenumber k% = k /a(T ) = 0.002 Mpc -1 [11], respectively. The 
corresponding physical frequency is /o — 3.09 x 10~ 18 Hz. For tensor power spectrum, one has the 
definition: 

Al(k) = h 2 (k,T ). (35) 

Hence, the non-zero value of r implies the existence of gravitational wave background, and would be 
probed with measurements of the B-mode CMB polarization [53[ ■ Even though there is still no direct 
observation of A^&nh A|(fc ) can be fixed by CMB observations, A^(k ) = A s = (2.43±0.11) ■ 10~ 9 
by WMAP 7 Mean |34lj . On the other hand, at present only observational constraints on r have been 



given by WMAP |33|, |34|. The upper bounds of r are recently constrained [34] as r < 0.24 by 



WMAP+BAO+ifo and r < 0.36 by WMAP 7 only for the vanishing scalar running spectral index 
a s , and r < 0.49 for the non- vanishing a s by both the combination of WMAP+BAO+ifo and the 
WMAP 7 only, respectively. Furthermore, using a discrete, model-independent measure of the degree 
of fine-tuning required, if 0.95 < n s < 0.98, in accord with current measurements, the tensor-to-ratio 
satisfies r > 10~ 2 jH3|. Therefore, one can normalize the RGWs at k = k using Eq. (13"4"|) . if r can 
be determined definitely. 

Since ku < &o < &2> h(ko,To) re-entered the horizon during the matter domination and is now 
inside the present horizon. Therefore, h(ko,To) suffered a decay oc 1/a from the horizon re-entry to 
the present time. According to Eq. (125]) one has 

h(k ,r ) = A h (k ) = A ( (1 + ze)-^. (36) 



H 



Therefore, A can be fixed for a determined value of A/ l (/c ), given a definite value of /3. Note that, the 
"thin-horizon" approximation that treats horizon re-entry as a "sudden" or instantaneous event 52| 
has been used in Eq. (l36l) . Thus, Eq. (l36|) is a normalization of RGWs provided by the observations of 
CMB. On the other hand, from the point of view of theories, the initial condition of RGWs could also 
be given, or some indirect connections between different parameters would exist. In the following, 
we discuss two different theories. 



Firstly, the constant A appearing in Eq. (|2~T|) can be determined by quantum normalization |2j: 

A = 68v^^, (37) 
where b = 7 2+/3 /|l + in our notation, Ipi = a/G is the Plank length, and Iq has been fixed to 



be 42 



/ = fel^'Cr* Cf C2 2 C E 7 , (38) 

using the continuous conditions of a(r) and a'(r) at the joining moments between two different 
phases. Here, it should be point out that the values of ( s and £i shown in Eq. (|38|) will be chosen 
tentatively as those determined by Eqs. (Tl0|) and ( ITT]) , respectively, since they are not known well. 
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FIG. 3: The relation between r and /3. Left: Three r — /3 curves based on the quantum normalization for three different values 
of the scalar spectral index n a — 0.967, n a — 0.966 and n s — 0.963, respectively. Right: The comparison between the r — ft 
relations based on the quantum normalization and those based on the slow-roll inflation. 



With the help of Eqs. (l3l), (l37j) and (l38|), Eq. (l36l) reduces to 

An(k y/ 2 = S^lpiHoCp^C'C^C^ (j^Y, (39) 

which implies a relationship between (3, /3 S and r. Given /3 a — 1, we illustrate r as a function of /3 
in the left panel of Fig. [3] for n s = 0.967, 0.966 and 0.963, respectively. Here we plus the case of 
n s = 0.963 for comparison. One can see that the r — (3 curves are almost overlap for (3 ^ —2 and 
only small discrepancies exist for lower (3. Taking n s = 0.966 for example, the constraints r < 0.49, 
r < 0.36 and r < 0.24 lead to (3 > -2.038, (3 > -2.035, and (3 > -2.032, respectively, while r > 0.01 
gives (3 ^ —2.005. Note that, these results are based on the validity of quantum normalization Eq. 
f )37|) . however, it is not the unique initial condition. This is a manifestation of the vacuum ambiguity 



that is responsible for particle production in cosmological spacetimes [55j. Furthermore, we have 
taken the values of ( s and £i evaluated by using the slow-roll inflation model with a scalar potential 
V(4>) = \w?^) 2 {4o[. If other different values of ( s and Ci are chosen, we find that the curves plotted 
in Fig. |3] will change sensitively. For simplicity, in this paper we will not consider these cases. 

Secondly, from the point of view of the slow-roll inflation model, there is a natural relation between 
n s and r. For V{4>) = ^m 2 (p 2 , using Eq. (jSJ) one easily has e = r/ and 

n,«l- 4e. (40) 

On the other hand, under the slow-roll approximation, the primordial tensor power spectrum and 
the primordial scalar power spectrum are respectively given as 



Al(k,n)*™(^-) 2 , (41) 
A 2 n (k,n)^-(^-) 2 , (42) 



where if* is the Hubble rate during the inflation stage and is invariable for the slow-roll approxi- 
mation, r* stands for the moment when the k-mode exits the horizon. According to the original 
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definition of r, i.e., the ratio of the primordial tensor power spectrum to the primordial scalar power 
spectrum, one has 45|, |52 



r = A 2 h (k,T*)/A* l (k,T*) = 16e, (43) 
with the help of Eqs. (I4ip and (142)) . Note that the definition of r in Eg. (I43p is little different from 



that in Eq. ( 154)) . We will discuss the difference below. In the literatures of WMAP |31|, |32|, [34|, the 
primordial power spectrum are often written as 

Al(k) = Al(k ) (j-J , (44) 
/ k \ ns_1 

AUk) = AUh) - , (45) 



where n t is the tensorial spectral index. Generally, n t is k— depedent [49l . |51| . l56l . |57J, however, for 
simplicity we just consider n t as a constant in this paper. With the help of Eq. ( 12T|) . one easily find 
that n t = 2/3 + 4. Strictly speaking, based on Eq. ( 126|) . Eq. (jUJ) is valid approximately since the 
mode k = k Q has re-entered the horizon at the present time. An analogous case would also exist in 



Eq. ( )45|) . which was used to evaluate the reheating temperature in [40|. However, as a conservative 
consideration, Eqs. ( 14*4)) and ( 145]) provide a good approximation to evaluate r. Now let us estimate 
the discrepancy between Eq. (IMj) and ( 1431) in the following. From Eq. (j4*4"j) . it is straight forward to 

S et r = ^iSy = ^§(W^o) nt+1 ~ ns - The differential factor (k E /k ) nt+1 - n ° ~ 0.992 for n s = 0.966 
or n s = 0.967, where we have used Eq. ( 140|) and the well known relation 

n t « -2e (46) 

in the slow-roll inflation. Allowing for the relation n t = 2/3 + 4, the combination of Eqs. (j4"3"j) and 
(1461) give 



r = -8n t = -1603 + 2). (47) 

This general linear relation is plotted in the right panel of Fig. |3j For a fixed value of n s , r can 
be determined through Eqs. ( 140]) and (143]) . and then /3 can also be determined through Eq. ( 147)) 
correspondingly. We employ a red circle and a green square to stand for the values of (/3, r) = 
(—2.009,0.136) and (/3,r) = (—2.008,0.132) corresponding to the observations n s = 0.966 and 
n s = 0.967, respectively. For comparison, the relation between r and (3 determined by the quantum 
normalization was shown again, and the constraint of 0.01 < r < 0.24 discussed above leads to /3 
localizing in the range of —2.032 < /3 < —2.005. On the other hand, according to Eq. ( 14 7p . the 
constraint of 0.01 < r < 0.24 forces /3 to lie in the range of —2.015 < /3 < —2.001. Hence, under the 
condition r > 0.01, /3 has a tilt smaller than —2 both for the two theoretical bases. In fact, Eq. ( 14T)) 
directly implies /3 < —2 due to the definition of r being positive. 

In the following, we demonstrate the properties of the energy spectra of RGWs taking some definite 
values of r and /3 for example. Firstly, setting n s = 0.966, the approximate energy spectra with three 
different tensor-to-scalar ratios r = 0.49, r = 0.01 and r = 0.001, which correspond to (3 = —2.038, 
f3 = —2.005 and (3 = —1.985 according to quantum normalization, respectively, are illustrated in 
the left panel of Fig. UJ It can be found that all the three cases of energy spectrum have the same 
amplitudes at / = /i ~ 8.04 • 10 7 Hz. This is a natural result from the quantum normalization. 
Plugging Eq. (ETJ into Eq. OH with the help of Eqs. ([23]) and ([38]), one obtain 

hc(fi) = tVtofrlpu (48) 
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and it is straight forward to obtain 

a a 3 fAl2 

n,(/0 = (49) 

according to Eq. (1321 . We can see that h c (fi) and fl g (fi) are fixed values independent on 0, 
since /i is independent on /3. On the other hand, as shown in Fig. HI a smaller (3 leads to larger 
energy spectrum especially at lower frequencies. This property is contrary to that illustrated in 



Refs. |41|, |42|, |51j, where a fixed (3 S was chosen but the quantum normalization was not employed. 
Secondly, let us see the differences between the energy spectra of RGWs based on the quantum 
normalization and those based on the slow-roll inflation. In the right panel of Fig. HJ we plot 
the energy spectra of RGWs based on the two different theories for n s = 0.966 and n s = 0.967, 
respectively. For the slow-roll inflation, once n s is given, r and (3 will be determined subsequently. 
Concretely, as analyzed above, n s = 0.966 and n s = 0.967 correspond to r = 0.136 and r = 0.132, 
respectively. For comparison, we also choose the same values of r in the energy spectra based on 
quantum normalization, and moreover, it is easy to have (3 ~ —2.027 for r = 0.136 or r = 0.132 
using Eq. ( 139]) . As shown in the right panel of Fig. [4] the energy spectra of RGWs based on these 
two theories are almost the same at / < f 2 both for n s = 0.966 and n s = 0.967. This is because of 
the normalization from the CMB observations at / = fo, where Q g (f)/r is a fixed value. Therefore, 
around / = f , fi g (f) depends on r sensitively rather than (3, since Q g (f) depends on j3 in the way 
of Sl(f) oc {f/U) w . Let us take fi 9 (/ 2 ) for example. With the help of Eqs. (EHD, fl32]) and ([36]), one 
can easily derive that VL g (f2) oc r (^) 2/3 ("^) 2 - Among the four curves illustrated in the right panel of 
FigJU the biggest relative discrepancy of Q g (f 2 ), existing between the sets (r = 0.136, (3 = —2.008) 
and (r = 0.132, (3 = —2.027), is about 17%. Similarly, for the modes which re-entered the horizon 
during radiation domination, one has Q g (f) oc r (^) 2/3 (^) 4 due to Eq. (T29]) . Therefore, for fixed r, 
the discrepancies of the energy spectra based on quantum normalization and those based on slow- 
roll inflation are larger and larger with the increasing frequency, since a difference A/3 ~ 0.02 in the 
indices between them always exists. In the very high frequency band f s < f < fi, not only the 
discrepancy discussed above is larger, but also the energy spectra manifest different properties for 
different values of n s . The latter phenomenon results from that n s = 0.966 and n s = 0.967 lead to 
different values of f s and fi as analyzed below Eq. (133]) . either on the basis of quantum normalization 
or slow-roll inflation. The value of f s depends on n s very sensitively, while f\ depends on n s much 
less sensitively compared to f s . These properties would become an effective tool to discriminate the 
value of n s , and could be an interesting object of the high-frequency gravitational wave detectors 

23411. 



Below, let us discuss the detection of RGWs using the ongoing and planed gravitational detectors 
which are sensitive at different frequency bands. As shown in the right panel of Fig. H] there is nearly 
no difference between the RGWs with n s = 0.966 and those with n s = 0.967 at / < f s , both based on 
quantum normalization and slow-roll inflation. Thus, we just take n s = 0.966 for demonstration in 
the following. As a conservative evaluation, in Fig. |5]we show the strain sensitivity curves of various 
gravitational wave detectors and the strain amplitudes, h c (f) / yff, of the RGWs for n s = 0.966 based 
on quantum normalization and slow-roll inflation, respectively. These detectors contain the PPTA 



25l | and SKA [58 J u sing the pulsar timing technique, th e sp ace-based laser interferometers such as 



LISA [20], BBO |21|,|22j], and the Fabry-Perot DECIGO [23j, the ground-based laser interferometers 
including the first generation LIGO [§| and VIRGO [l3j], the second generation AdvLIGO 11] and 



LCGT [16j, and the third generation ET [18|. One can see that all the ground-based interferometers 
can hardly detect the theoretical RGWs discussed above. PPTA and LISA also have difficulties to 
catch RGWs, however, the planed SKA, BBO and DECIGO are promising to detect RGWs since 
they have relative higher sensitivities. In order to show the detectabilities of RGWs by SKA and 
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FIG. 4: Left: The energy density spectra of RGWs based on the quantum normalization with r = 0.49(/3 = —2.038), r = 
0.01(/3 = —2.005) and r = 0.001(/3 = —1.985), respectively, for n s = 0.966. Right: The comparisons between the energy density 
spectra based on the quantum normalization and the slow-roll inflation for n s — 0.966 and n s — 0.967, respectively. 




10-" 10- 5 10° 10 5 

Frequency [Hz] 



FIG. 5: The strain of RGWs with different parameters for n 3 = 0.966 confronts against the current and planed GW detectors. 
The sensitivity curves of PPTA and SKA using pulsar timing technique are taken from Refs.[25| and [58|. respectively. The 
curve of BBO is generated using the online "Sensitivity curve gene rator" [2p| ] with the parameters in Table II of Ref. [2l[ and 
Table I of Ref. [22]. The curve of FP-DECIGO is taken from R efTpljl . The sensitivity curves of the ground-based interferometers 
AdvLIGO, LCGT and ET are taken from Refs.QjJl, Q]| and [irrespectively. 



DECIGO/BBO more clearly, we enlarge the two parts of SKA and DECIGO/BBO which are shown 
in Fig. [HI As can be seen in the left panel of Fig. El around the frequencies 3 — 5 x 10 -9 Hz, the 
spectra of RGWs with different parameters of r (and in turn the different corresponding /3 due to the 
r — relation) based on quantum normalization have different amplitudes, and they are promising to 
be distinguished by SKA. However, the spectrum based on slow-roll inflation is hardly distinguished 
from that based on quantum normalization with r = 0.49(/3 = —2.038) by SKA, since two lines 
almost overlap in the frequency response band of SKA. On the other hand, in the right panel of 
Fig. [6] one can see that, RGWs based on slow-roll inflation are easier to be distinguished from those 
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FIG. 6: A zoom plot of the two parts of SKA and FP-DCIGO/BBO. All the values of parameters in RGW models and the 
sensitivity curves of the three gravitational wave detectors are exactly the same as those shown in Fig. [5] 



based on quantum normalization by DECIGO or BBO in the frequency range ~ 10~ 2 — 10° Hz. 
Furthermore, since DECIGO/BBO has a wider frequency response band in contrast to SKA, the 
r — (3 has the potential to be examined. In conclusion, the combination of SKA and DECIGO/BBO 
provides an important tool not only distinguishing different theoretical models of RGWs but also 
determining the corresponding parameters. 



V. DISCUSSIONS 



We determined the expansion histories of the preheating stage and the radiation-dominated stage 
due to the the reheating temperature which is recalculated using the latest observations of WAMP 
7- year. Based on that, we illustrated the approximate solutions of RGWs in the current accelerating 
stage for various frequency bands. We found that the frequency f s , describing that the mode re- 
entered the horizon at the moment of the reheating, is dependent on n s sensitively, however, the 
upper limit frequency of RGWs depends on the value of n s much less sensitively. Combing the 
quantum normalization of RGWs with the CMB observations, we obtained a relation between the 
tensor-to-scalar ratio r and the inflation index (3 for the fixed preheating index (3 S = 1. According 
to the relation between r and (3 with a fixed n s = 0.966, we find that a relatively tight constraint 
0.01 < r < 0.24 leads to (3 localizing in a range of —2.032 < (3 < —2.005 based on quantum 
normalization and in a more narrow range of —2.015 < (3 < —2.001 based on slow-roll inflation, 
respectively. We plotted the spectrum and the energy density spectrum of the RGWs for three cases 
of r = 0.49(/3 = -2.038), r = 0.01(/3 = -2.005), and r = 0.001(/3 = -1.985), respectively. It 
was found that a lager r, i.e., a smaller (3 leads to a larger spectrum of RGWs especially at lower 
frequencies. For comparison, we also illustrated the spectra of RGWs with the parameters of r and 
f3 given by the slow-roll inflation. Concretely, for n s = 0.966, one has r = 0.136 and [3 = —2.009, 
and for n s = 0.967, one has r = 0.132 and (3 = —2.008. It was found that, for the same value of r, 
the discrepancy of the energy spectra based on quantum normalization and those based on slow-roll 
inflation is larger and larger with the increasing frequency. However, our analysis above does not not 
apply, in general, to less conventional models of inflation where the RGW spectrum and the observed 
spectrum of scalar perturbations are produced with different primordial mechanisms 59, 6(|, and, 
as a consequence, they are in principle completely decoupled. 
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Among the current and planed GW detectors, only the planed SKA using the pulsar timing 
technique, and the planed space-based interferometers BBO and DECIGO are promising to detect 
RGWs. However, these results are based on the referenced values of £i and £ s which are obtained 
from the combination of CMB observations and the slow-roll inflation with a concrete potential 
V{4>) = ^m 2 (j) 2 . Note that, the values of (j and ( s depend sensitively on the reheating temperature 
Trh which is dependent on n s sensitively. However, the r — (3 relation is nearly not dependent on n s 
once the form of the potential V(<f>) = |m 2 2 is given. Hence, for the same value of r, different n s 
only affects the very high frequency RGWs that re-entered the horizon before the radiation-dominant 
stage. Therefore, to a certain extant, our analysis of the detection in the frequency range / < 10 3 
Hz is general, even though the determination of the reheating temperature from CMB has very large 
uncertainties. As forecasted in Ref. [401. the future CMB experiments such as the Plank satellite 



35[, the ground-based ACTPol |36| and the planned CMBPol [37( will provide significant reductions 
of the uncertainties of the reheating temperature Trh- Therefore, one expect accordingly that the 
expansion histories of the very early universe would be known better. On the other hand, the values 
of £i and ( s could be chosen independently on the slow-roll inflation scenario, which would be studied 
in the future work. 
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